First principles calculation of structural and magnetic properties for Fe monolayers 
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Structure optimizations were performed for 1 and 2 monolayers (ML) of Fe on a 5 ML W(llO) 
substrate employing the all-electron full-potential linearized augmented plane-wave (FP-LAPW) 
method. The magnetic moments were also obtained for the converged and optimized structures. 
We find significant contractions (~ 10 %) for both the Fe-W and the neighboring Fe-Fe interlayer 
spacings compared to the corresponding bulk W-W and Fe-Fe interlayer spacings. Compared to the 
Fe bcc bulk moment of 2.2 /xs, the magnetic moment for the surface layer of Fe is enhanced (i) by 
15% to 2.54 /xs for 1 ML Fe/5 ML W(llO), and (ii) by 29% to 2.84 /xs for 2 ML Fe/5 ML W(llO). 
The inner Fe layer for 2 ML Fe/5 ML W(llO) has a bulk-like moment of 2.3 ^b- These results agree 
well with previous experimental data. 
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I. INTRODUCTION 

Magnetic thin films on metal substrates demonstrate 
fascinating phenomena such as the preferential orienta- 
tion of the magnetization normal to the film plane, en- 
hanced low-temperature surface magnetization, and the 
pronounced effects of magnetism on the electeical conduc- 
tivity. There has been considerable efforttrEll in studying 
the atomic structures and magnetic properties including 
magnetic moments and the orientation of the magnetic 
easy axis for Fe thin films on W(llO). It is especially 
interesting to study these properties for 1 and 2 mono- 
layers (ML) of Fe on W(llO) due to the pseudomorphic 
layer-by-layer growth of the film when the Fe coverage 
9 is below 2 ML and to the possible magnetization re- 
orientation from in-plane to perpendicular for 1 < 6 < 
2. Both bulk Fe and W are bcc structures with lattice 
constants of 2.86 A and 3.165 A respectively. It was 
found that-Fe thin films grow pseudomorphically up ti 
1.2 MLBEI on the flat W(llO) surface and up to 1.8 MlI 
on a vicinal surface. Significant structural relaxation in 
the vertical lattice spacings for the Fe thin films is ex- 
pected arising from the large lattice mismatch (~ 9 %) 
between the film and the W substrate. However, the ex- 
act amount of relaxation and the magnetic moments for 
the Fe overlayers need to be clarified because-pf the con- 
flicting results between the expjerimentsEra'Ej and with 
previous theoretical predictionsu. On account of (i) the 
pseudomorphic growth of the Fe thin film, and (ii) the 
transitions of both the atomic structures and magnetic 
properties already in the ultrathin (<2 ML) regime, it 
is feasible to employ the ab initio method to investigate 
these properties. __ 

Earlier experimental work done by Albrecht et aluu 
with low-energy electron diflraction (LEED) on 1 ML 
Fe on W(llO) substrate showed that the Fe-W inter- 
layer spacing is contracted by 13% to 1.94 A compared 
to the bulk W(llO) interlayer spacing of 2.238 A. The 



magnetic moment for the top Fe layer is enhanced to 
2.53 /ifl-jjncasured by Torsion Oscillation Magnetometoj 
(T0M)Ei3. However, recent work done by Tober et a£j 
using Photo-Electron Diffraction (PED) for 1 ML Fe on 
W(llO) yielded a Fe-W interlayer spacing of 2.07 A, a 
relaxation a£ 7.2% only. Earlier ah initio calculations by 
Hong et ala showed that the Fe-W interlayer distance 
is dramatically reduced by as much as 16% to 1.88 A. 
Because of this strong inward relaxation, the magnetic 
moment of the over layer Fe is only 2.18 /is which is very 
close to the bulk bcc Fe moment of 2.2 /i^. Recent cal- 
culations by Batirev et a/.E3 showed that the Fe-W inter- 
layer spacing is contracted by 3.1% with respect to the 
average theoretical bulk bcc(llO) Fe and W interlayer 
spacings. The magnetic moment for the Fe layer is 2.17 
/is similar to the bulk value. Because of these differences 
between the experimental results and also with the calcu- 
lations, a detailed theoretical investigation is presented 
to illuminate the incoherent data on structural and mag- 
netic properties of Fe monolayers on W(llO). 



II. METHOD 

It is well known that magnetic properties depend 
strongly upon the atomic structures of the thin films. 
Therefore it is necessary first to obtain an optimized 
structure for these systems. Three different slabs: (i) 
5 ML W(llO) clean substrate, (ii) 1 ML pseudomorphic 
Fe overlayer on each side of 5 ML W(llO) substrate, and 
(iii) 2 ML pseudomorphic Fe overlayers on each side of 5 
ML W(llO) substrate were studied. The schematic pic- 
ture of 2 ML Fe on 5 ML W(llO) is shown in Fig. 1. The 
bare W(llO) substrate was studied to test our theoretical 
accuracy since reliable experimental results are available 
and theoretical calculations are abundant. 

These calcvilations were performed employing the 
WIEN97 codecJ. This program is based on the density- 



functional-theory (DFT) and adopts the fuU-potential 
Unearized augmented plane-wave (FP-LAPW) method. 
It has-|the additional capability of computing atomic 
forcescScS which makes the structure optimization much 
more efficient compared to the total energy only calcu- 
lations. The FP-LAPW method adopts different repre- 
sentations for wave functions, charge density and poten- 
tial inside the muffin-tin sphere and in the interstitial 
region. The spherical harmonics were expanded up to 
1=10 inside the muffin-tin and to 1=4 for the intersti- 
tial in the present calculations. Spin-polarized calcula- 
tions were carried out in order to determine the magnetic 
properties. Spin-polarization was implemented in the 
WIEN97 code adopting Local Spin-Density Approxima- 
tion (LSD A) with two separate spin densities. Two sets 
of Kohn-Sham (KS) orbitals for the two spin components 
were obtained, and two sets of KS single particle equa- 
tions were solved. The scalar relativistic calculations in- 
cluding the velocity and the Darwin terms were adopted 
for valence electrons. Spin-orbit coupling for the valence 
electrons was not included for the present calculations. 
The shallow 5p states were treated as semi-core, i.e. as 
local orbitals, thereby ensuring the flexibility of the basis 
functions to closely^i^epresent these low-lying p orbitals. 
Pulay correctionsE^EZI to the Hellmann-Feynman forces 
were calculated which makes the structure optimizatiou, 
highly accurate. The improved tetrahedron methodES 
was used for the integrations. 

Fig. 2 shows the unit cells for the calculations. The 
structure optimization for the slabs was done by giving 
an initial guess of the interlayer spacings based on the 
optimized structure of Fe/Mo(110). The direction and 
degree of relaxation for the vertical interlayer spacings 
depend on the magnitude and sign of the forces present. 
The in-plane lattice constant for the slab was fixed and 
taken from the bulk calculations and will be described 
later. This is due to the in-plane two-dimensional trans- 
lational invariance and the fact that there is only one 
atom on each layer in the unit cell. As described in our 
previous papeicJ, eight vacuum layers were incorporated 
in the supercell to separate the slabs in order to mini- 
mize any Coulomb and exchange interactions. Further- 
more, slabs are symmetric with respect to the central sub- 
strate layer to avoid any charge accumulation on the sur- 
faces. Thus the contribution to the total energy from the 
electric-dipole interaction between the supercells is neg- 
ligible compared to the contributions from within the su- 
percell. In addition, only real wave functions are needed 
for the calculations because of the presence of inversion 
symmetry. The Fe layers on each surface are ferromag- 
netically coupled. The spin-polarized calculations were 
applied. 

In these calculations, Generalized Gradient Ap- 
proximation (GGA)E3 exchange potential and scalar- 
relativistic treatment were used in agreement with our 
earlier calculations on Mo substrate. Generally speak- 
ing, we did not find any significant improvement of GGA 
exchange potential over LSDA potential. Following the 



procedure described previouslyE^, the theoretical bulk 
W lattice constant was determined to be 3.205 A, 1.3% 
larger than the experimental value of 3.165 A. It is known 
that GGA corrects overbinding, but sometimes leads to 
an excessive increase in the lattice parameter for heavy 
atoms such as W. Nevertheless this theoretical value was 
used as the in-plane lattice spacing in our subsequent 
slab calculations. The theoretical bulk Fe lattice con- 
stant was found to be 2.834 A, 0.9% smaller than the 
experimental result of 2.86 A. The muffin-tin radii were 
chosen to be 1.27 A and 1.164 A for W and Fe atoms 
respectively in the slab unit cells. Convergence was 
achieved when the total energy and charge differences 
between two consecutive iterations are less than 5xl0~^ 
Ry and IxlO"** e/(a.u.)'^ respectively. The structure op- 
timizations were done when the force on each atom is 
less than 1 mRy/a.u.. The magnetic moments were cal- 
culated as the differences between the spin-up charge 
and spin-down charge for these converged results. Or- 
bital magnetic moment is not included in our calculations 
due to the absence of spin-orbit coupling for the valence 
electrons. Moipepver, it was previously estimated to be 
around 0.1 /is" only. The numbers of k points in the 
two-dimensional meshes are 20 x 20 for 5 ML W(llO), 
21 X 21 for 1 ML Fe on 5 ML W(llO), and 22 x 22 for 
2 ML Fe on 5 ML W(llO). The numbers of k points 
in the irreducible part of the Brillouin zone (IBZ) (1/4 
of BZ) are 110, 121, and 132 respectively. The plane- 
wave cut-offs (corresponding to the largest k-vector in 
the plane- wave basis expansion) are 16.7, 15.3, and 13.2 
Rys for the three slabs respectively with 0, 1 and 2 ML 
Fe coverage. The kinetic energy cut-offs (corresponding 
to the largest reciprocal-space vector for the potential 
expansion) are 196 Ry for all three slabs. 



III. RESULTS AND DISCUSSION 



The structural and magnetic results are exhibited in 
Tables I and II respectively. For the 5 ML W(llO) clean 
substrate, we find that the top W-W interlayer spacing 
is contracted by 4.1% to 2.173 A from the theoretical 
bulk W-W interlayer spacing of 2.266 A in the (110) 
plane. This result is in sopd agreement with the previous 
FP-LAPW calculations^ in which the same amount of 
contraction (4.1%) was found for the top W-W interlayer 
spacing with a 5 ML W(llO) slab. A 3.6% downward re- 
laxation was found for the top layer with a 9 ML W(llO) 
slab. iQur result is in disagreement with a recent calcu- 
lationca. However, in that study, only three substrate 
W(lip)| layers were employed. The recent LEED exper- 
iments yielded a contraction of 3.1% with an error bar 
of 0.6%. In addition to our agreement with previous the- 
oretical and experimental data, our present result is also 
quite similar to the redaxation found for a 5 ML Mo(llO) 
slab published earlieiE3. Further, we find that the sec- 
ond W-W interlayer distance is also slightly contracted 
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by 0.4% tO|2,258 A. This again agrees well with earlier 
calculationsEj in which a 0.2% contraction was found for 
the 5 ML W(llO) slab. 

The clean W(llO) substrate is non-magnetic. The 
density-of-state (DOS) plot is shown in Fig. 3. Only d- 
partial DOS (PDOS) of spin-down are shown since they 
are identical to the spin-up DOS. The inner W layer 
(W(S-2)) rf-PDOS closely resembles the bulk bcc W d- 
PDOS. The surface layer (W(S)) rf-PDOS has a higher 
number of states at the Fermi-level, almost double that 
of the W(S-2), i.e. a less pronounced gap between the 
two sub-bands. 

For the slab of 1 ML Fe on each side of 5 ML W(llO) 
substrate, we find a significant relaxation for the Fe-W 
interlayer spacing (Table I) very similar to the case otJ 
ML Fe/5 ML Mo(llO) as shown in our previous worked. 
The Fe-W interlayer has a downward relaxation of 12.9% 
compared to the bulk W-W interlayer distance. It is in 
excellent agijeepent with the LEED experiment by Al- 
brecht et aluU in which a 13% contraction was found 
compared to the bulk experimental W-W interlayer dis- 
tance. The recent PED experiment, however, yielded a 
Fe-W distance of 2.07 A(7.2% contraction only) which 
corresponds to the bond length from the hard sphere 
model. The earlier calculations done by Hong et al 
showed a much larger downward relaxation of 16% em- 
ploying the FP-LAPW method. However, in their earlier 
calculations, not all the atoms were allowed to relax at 
the same time since it was not possible to compute the 
force on each atom. In addition to the Fe-W distance, 
our present calculations show that the neighboring W-W 
interlayer spacing is reduced slightly by 0.1%. However 
our earlier results on Mo(llO) show a small expansion 
for the neighboring Mo-Mo interlayer spacing contrary 
to the W case here. 

The magnetic moment for the surface layer of Fe is 
found to be 2.54 ^b without orbital moment contribu- 
tion, an enhancement of 15% over the bulk magnetic mo- 
ment of 2.2 jiB for bcc Fe. However, it is reduced by 29% 
compared to the moment of 3.3 ^b for the Fe(llO) free- 
standing monolayer with the same in-plane lattice pa- 
rameter. In addition, our results show that the neighbor- 
ing W layer acquires a small moment of 0.1 ^b- It is anti- 
ferromagnetically coupled to the Fe overlayer. The Tor- 
sion Oscillation Magnetometcy(TOM) experiment done 
by Gradmann and coworkersli3 yielded a moment of 2.53 
/is for the overlayer Fe. Since the orbital moment and 
the induced substrate moment are both around 0.1 /is 
and opposite in sign, the theoretical spin moment we ob- 
tained for Fe overlayer agrees very well with the TOM 
experiment since TOM measures the total mom£nt. The 
earlier calculations done by Hong and coworkerto showed 
no enhancement of the moment over the bulk value. It 
is probably due to the fact that their calculations yield a 
significant reduction of the Fe-W interlayer distance. 

Spin-down and spin-up d-PDOS for both the surface 
Fe layer and the neighboring W layers are plotted in Figs. 
4 and 5 respectively. The d-PDOS of both spins for the 



inner W(S-2) layer are very similar to the central W layer 
of the previous case where there is no Fe overlayer. The 
spin-down rf-PDOS of the interfacial W(S-l) layer resem- 
bles the one of the inner W(S-2) layer. However, there are 
noticeable changes for the spin-up component especially 
when close to the Fermi surface. The small moment of 
the W(S-l) layer is due to this change of d-PDOS. The 
overlayer Fe d-PDOS are very different to their corre- 
sponding bulk ones as shown in Figs. 6 and 7, especially 
for the spin-up component. Figs. 6 and 7 will be dis- 
cussed a little later. 

For the slab of 2 ML Fe on each side of 5 ML W(llO) 
substrate, we find both the Fe-Fe and Fe-W interlayer 
spacings are contracted dramatically (see TaMe I). The 
Fe-Fe interlayer distance is reduced by 11.9 %E^ from the 
theoretical bulk Fe value of 2.004 A to 1.766 A. The Fe- 
W interlayer spacing is contracted by 10.6% compared to 
the bulk W-W interlayer distance. The percentage of the 
contxactions are also vgry close to the Fe(110)/Mo(110) 
casec^. Albrecht et a/.Era found a 10% downward relax- 
ation for the Fe-Fe interlayer spacing compared to the 
bulk Fe value. Again it is in excellent agreement with 
our findings. Our calculations show a slight expansion < 
0.2% for the inner W-W interlayer spacings. 

The magnetic moment for the surface layer of Fe is 
found to be 2.84 /is, an enhancement of 29% over the 
bulk value of 2.2 jjlb- It is still smaller than the moment 
of the Fe(llO) free-standing monolayer. However, com- 
pared to the 1 ML Fe/5 ML W(llO) case, the moment 
for the top Fe layer is increased from 2.54 ^.b to 2.84 
/is. This is probably due to the strong hybridization of 
Fe d-orbitals with the ones of interfacial W layer thereby 
reducing the moment of the Fe layer. The second Fe 
layer, i.e. the interfacial Fe layer has a moment of 2.3 
/is already very close to the bulk value. Like the pre- 
vious case, the neighboring substrate layer also acquires 
a small moment of 0.1 hb and is antiferromagnetically 
coupled to the Fe over layers. 

The d-PDOS for the surface and interfacial Fe layers 
together with the ones of Fe bcc bulk are plotted in Fig. 
6 (spin-down) and Fig. 7 (spin-up) for comparison. Ba- 
sically the d-PDOS of the second layer of Fe are already 
close to the bulk ones. Consequently its magnetic mo- 
ment is also approaching the bulk value. The d-PDOS 
of the surface Fe layer are different from the bulk ones 
particularly for the spin-up component and when close 
to the Fermi-surface for the spin-down component. 



IV. SUMMARY 

The present FP-LAPW calculations resolves the dis- 
crepancies between previous experimental data and with 
earlier theoretical results on the atomic structure and 
magnetic moment of 1 ML Fe/W(110). The Fe-W in- 
terlayer spacing is significantly contracted by as much 
as ~13% compared to the bulk W-W interlayer spacing. 



The magnetic moment of the overlayer Fe is greatly en- 
hanced compared to the bulk moment of bee Fe due to 
the lower coordination number, but it is reduced com- 
pared to the Fe(llO) free-standing monolayer because of 
the presence of the substrate. 
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d(Fe2--Fei) 


d{W2-Fei) 


d{Wl-W2) 


d(W2-W3) 


5 ML W(llO) 

1 ML Fe/W(110) 

2 ML Fe/W(llOi 


L766(-11.9%)EI 


L974(-12.9%) 
2.026(-10.6%) 


2.173(-4.1%) 
2.263(-0.1%) 
2.267 (0.03%) 


2.258(-0.4%) 
2.251(-0.7%) 
2.272(0.2%) 


W(llO) (Exp.)yp 
Fe/W(110)(Exp.)2 
Fe/W(110)(Exp.)i, 
Fe/W(110)(Exp. ia 
1 ML Fc/W(110)'y 


L82(-10%) 


L94(-13%) 
2.07(-7.2%) 
L88(-16%) 


2.169(-3.1%) 
2.28(2.2%) 





TABLE L Structural results (The layer spacings are given in A. The relative changes as compared to the bulk W layer 
spacing are given in parentheses. The percentage of Fe-Fe contraction is relative to the bulk Fe-Fe interlayer spacing.) 





1 ML Fe/W(110) (/is) 


2 ML Fe/W(110) (/is) 


Fe(2) 




2.850(2.844) 


Fe(l) 


2.561(2.536) 


2.315(2.308) 


W(l) 


-0.085 


-0.104 


W(2) 


-0.000 


-0.004 


W(3) 


-0.000 


-0.006 


Interstitial ^ 


-\ 


-0.055 


-0.04 


Fe(2)(Exp.)Fe/W(110))il 




2.77 


Fe(l)(Exp.)Fe/W(110))t3 


2.53 





TABLE II. Magnetic moments (values in the parentheses are obtained by adding interstitial contributions). 
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FIG. 1. Schematic picture of 2 ML Fe on 5 ML W(110)(upper half of the slab only). 
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FIG. 2. Unit cells of the 1 ML Fe on each side of 5 ML W(llO). 
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FIG. 3. Spin-down partial-rf density-of-states(d-PDOS) for 5 ML W(llO) clean substrate. S represents the surface layer, S-1 
the layer next to the surface layer, and S-2 the central layer. 
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FIG. 4. Spin-down partial-d density-of-states(d-PDOS) for 1 ML Fe/5 ML W(llO). 
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FIG. 5. Spin-up partial-d density-of-states(d-PDOS) for 1 ML Fe/5 ML W(llO). 
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FIG. 6. Comparison of the Fe spin-down partial-d density-of-states(ii-PDOS) for 2 ML Fc/5 ML W(llO) with the corre- 
sponding bulk bcc Fe one. 
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FIG. 7. Comparison of the Fe spin-up partial- lidensity-of-states(ii-PDOS) for 2 ML Fe/5 ML W(llO) with the corresponding 
bulk bee Fe one. 
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